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As oscillator applications, resistance-inductor-capacitor (RLC) circuits are 
employed in a diversity of settings. A low-pass, band-stop, band-pass, or 
high-pass filters can all be designed using an RLC circuit. A two-loop RLC 
circuit could not be represented mathematically in prior studies. Laplace 
transform is one type of integral transformation, which is able to resolve 
both second order non-uniform and uniform linear differential equations. 
This work solves the differential equations (DEs) of a two loops RLC circuit 
of an alternating voltage source by using two alternative approaches, 
Laplace transform (LT) and deep learning convolutional neural network 
(DLCNN). Initially, two DE have been declared. Next, Laplace transform is 
computed to solve these equations with symbolic variables for the first loop 
current and capacitor charge. Finally, we substitute the numerical values of 
the circuit elements for the symbolic variables. The charge and current 
initially decline exponentially. On the other hand, they oscillate over a long 
period of time. The capacitor charge and current initially decline 
exponentially and oscillate over a long period of time. The qualities of the 
result can be examined with a symbolic result, which is not possible with a 
numeric result. 
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1. INTRODUCTION 


Due to experimental or numerical analyses on the behavior between a dissipative and conservative 
system and the absence of fractionalized systematic techniques, recent research has been diverted from the 
importance of the recent fractional derivatives including the non-singular kernel with non-locality and the 
singular kernel with the locality. The mathematical representation of an ordinary differential equation (DE) 
with second-order cubic nonlinearity is the Van-der-Pol equation. The Van-der-Pol equation has been given a 
time delay in several investigations. Theresistor-inductor-capacitor (RLC) circuit differential equation is 
derived as a delay differential equation in this study together with the Van der Pol model differential 
equation [1]. Analytical solutions for the Caputo-Fabrizio, Liouville-Caputo, and new Mittag-Leffler 
function-based fractional derivative to describe the electrical RLC circuit model were previously discussed. 
The fractional differential equations take different sources into account. When the fractional order equals 1, 
the conventional behaviors are restored [2]. Dynamical system approaches and Melnikov theory can be used 
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to examine tiny amplitude perturbations of some implicit differential equations appearing in RLC circuits [3]. 
A steady-state process in an RLC circuit with power sources operating at unrelated frequencies is also 
considered [4]. An expansion of an ordinary differential equation is taken into consideration in order to 
achieve the periodic steady-state behavior. This expansion is based on changing from ordinary differential 
equations to partial differential equations with two-time variables by adding an additional time variable. The 
two-dimensional Laplace transform is used to solve the obtained differential equations. Active power and 
frequency responses for the domain of two time variables are specified by the use of double integral formulas 
for a transfer function. The voltage and current amplitude-frequency properties of the RLC circuit can be 
given in the domain of two variables. A nonlinear fractional derivative based Volterra integral-differential 
equation with Caputo, several kernels, and numerous constant delays is also considered to look into the 
qualitative properties of solutions to this equation, including the boundedness of nonzero solutions and the 
Mittag-Leffler stability, uniform stability, and asymptotic stability of the zero solution [5]. The 
Lyapunov-Razumikhin approach and selecting an acceptable Lyapunov function are the methods employed 
in the proofs of these theorems related to an RLC circuit. 

A modified Laplace transform approach to find solutions to a series-connected simple electric 
circuits (RLC) model of linear differential equations (DEs) was presented in [6]. Although the study 
suggested that non-homogenous second order linear differential equations in the form of electric charge 
equations can be solved over time using a modified Laplace transform method, only homogeneous second 
order linear DEs was considered. Certain approximation techniques for generated nonlinear terms of 
characteristic exponents to offer the stability analysis of delay integral-differential equations with fractional 
order derivatives were discussed in [7]. Such methods allow for the proof of the presence of some analytical 
solutions close to their equilibrium points. The formulas for the critical time delay and critical frequency are 
determined through the construction of stability charts to describe general RLC circuits that expose the delay 
and fractional order derivatives. A straightforward method was used to determine the electromagnetic energy 
density distribution in a dispersive and dissipative met material made up of wire arrays and split-ring 
resonators [8]. The paper demonstrated that the system’s energy may be thought of as being made up of the 
energy densities of the electric and magnetic fields as well as energy densities associated with the medium’s 
reaction. Therefore, the system’s equations of motion for polarization have been compared with the 
corresponding differential equations of suitable RLC circuits to develop formulas for the energy density of 
the medium. The electrical RLC circuit with a fractional DE is investigated in [9], where a new auxiliary 
parameter was added to maintain the three physical measures C, L, and R’s three dimensions. Through the 
circuit’s physical properties of RLC, the analytical solution was specified in terms of Mittag-Leffler’s 
formula. An analytical approach to examine the impacts of contemporary fractional differentiation on the 
RLC electrical circuit was discussed in [10]. The leading RLC electrical circuit DE has been fractionalized 
with three kinds of partial derivatives, where the RLC electrical circuit was seen for unit step sources, 
periodic, and exponential. 

Melnikov’s theory for implied ordinary DEs of little amplitude perturbations was presented in [11]. 
If certain Melnikov-like criteria are met, it is specifically explored how long orbits linking singularities can 
endure in finite time for nonlinear RLC circuit systems. The two delays in the Van der Pol delay model and 
delay differential equations were produced from ordinary differential equations using the Taylor series to 
define the model for treating Parkinson’s disease [1]. Floquet’s theory to examine the linear dynamic analysis 
of RLC circuits representing AC generators with periodically time-varying inductances was presented 
in [12]. The dynamic stability’s prerequisites are derived. The stability domains and transition curves were 
predicted using the harmonic balancing approach. A fractional order differential equation can be used to 
describe the dynamics of noninteger order RLC electrical circuits by adding an auxiliary parameter [13], 
where an analytical determination of the filter parameters and numerical validation of the outcomes were 
used. Despite the fact that differential equations are marketed as applied mathematics, the course rarely 
includes any practical applications [14]. Deep learning convolutional neural networks (DLCNN) [15]-[18] 
and reinforcement learning [19], [20] that were undergone fast growth in modern decades, were employed to 
solve the nonlinear DEs in general. Deep learning-based data assimilation algorithms [21], [22] were lately 
presented to train Navier-Stokes network formulas to estimate different quantities of interesting. 
DLCNN-based methods can be classified into three kinds: i) DLCNN that maps indirectly to the parameters 
or inner results of an algebraic explanation to use them for deriving numerical solutions [23], [24]; ii) 
DLCNN that maps directly toward the solutions with a discrete method, which is same as numerical solutions 
[25], [26]; and iii) DLCNN that maps straight to the solution characterized by a DLCNN with a continuous 
way and it is similar to that in analytical solutions [27]-[29]. In this type, the data applied for training the 
network are arbitrarily modeled inside the whole solution range in every training batch, including boundary 
and initial circumstances. The necessary characteristic of all of these techniques is the adopting the benefit of 
the nonlinearity depiction capability of DLCNNs. New development [30] in mechanics with this capability 
has been stated such as, Li et al. [31] presented a generative adversarial model network for mapping hidden 


Modeling two loops RLC circuit AC power source using symbolic arithmetic ... (Inaam Rikan Hassan) 


492 o ISSN:2302-9285 


variables for a microstructure to use in fabrication. The major contribution of this work is to model the DEs 
of a two loops RLC circuit of an alternating voltage AC-source by using two alternative approaches; Laplace 
transform and DLCNN in MATLAB environment. To obtain this aim, several objectives are accomplished 
such as: i) to model and solve the DEs of a 2-loop RLC AC circuit; ii) to plot and analyze the obtained results 
of the capacitor charging and inductor current; and iii) to solve and compare the results with the 
DLCNN-based solution. 


2. METHOD 
2.1. Object and research hypothesis 

This work solves the RLC circuit differential equations by using Laplace transform. It is possible to 
define the Laplace transform by a function f (t), which is given by: 


f roesat 
0 


Calculations are kept in their native symbolic form rather than in numerical form using symbolic 
workflows. This method enables using precise symbolic values and comprehends the properties of the 
obtained answer. When a quantitative result is required or symbolically continuing is impossible, numbers 
are substituted for symbolic variables. In general, the symbolic workflows of solving equations include; 
equations definition, equations solving, values substitution, plotting results, and analyzing results. 


2.2. Equations definition 

To solve differential equations using beginning conditions, utilize the Laplace transform. We will 
consider solving the RLC circuit shown in Figure 1. Capacitor charge is referred to as Q(t) in coulomb, AC 
voltage source in volts is referred to as E(t), capacitance in farad is referred to as C, inductance in henry is 
referred to as L, currents in ampere are referred as I1, I2, I3, resistances in ohm are referred as R1, R2, R3. 
Apply Kirchhoff s current and voltage laws, following equations, are obtained: 


h = l + Tz 
L% + hR, + 1,R, =0 (1) 
E(t) + hR, —2—1R3 = 0 


Figure 1. RLC circuit analysis 


By substituting the term 13 = dQ/dt (which represents the rate of charge of the capacitor) into 
eqn(1), the following differential equations are obtained for the RLC circuit: 


dl, R dQ R,+R, 


dt Ldt LoS 

(2) 
dQ g Ro 
dt Ro+R3 (EO z) t R2+R3 h 


Since the magnitude of the practical components are with positive quantities, we stat the variables of 
equations. Then, we set the equivalent statements for the variables. Suppose that E(t)is a 1 V AC source 
voltage (1. sin(t)), and assuming t, L, C, and R are greater than 0. Therefore, the differential equations are: 


Bulletin of Electr Eng & Inf, Vol. 13, No. 1, February 2024: 490-498 


Bulletin of Electr Eng & Inf ISSN:2302-9285 m) 493 


a 
a R2572 (t) 1 (t)(R1+R 
eqni(t) = 2h (t) - 2 = — Ae (3) 
P Q(t) 
a — Oey Rala(t) 
eqn2(t) = at Qg) = R2+R3 R2+R3 (4) 


2.3. Deep learning convolutional neural network framework 

From a physical perspective, the initial and boundary circumstances are continued with time, and as 
a result, the parameters of the approach circuit are also continuous in time: the parameters at the current step 
are used as an ideal initialization for the next step, allowing the parameters of the DLCNN to transmit with 
the time phase in the solution procedure, which makes the process rapid for the remaining time stages. The 
derivatives of the current and capacitor charge equations are computed according to the gradient of the 
network output according to the input. The hidden layers activation functions are selected as rectified linear 
unit (ReLU) functions to conserve the stability of those derivatives. The network architecture is trained 
according to the computed loss based on the gradient descent process. A common nonlinear differential 
equation has the following format: 


u, + N(u,9) =0 (5) 


The term u(x, t) represents the underlying equation solution, u; represents the time derivative, while 
N (u, 9) denotes the nonlinear function with parameter 0. Specifying the required condition, represented by 
the initial and boundary parameters, the equation solution is a nonlinear map within the solution range (x, t). 
DLCNNs have revealed noteworthy achievements in the learning of nonlinear high-dimensional 
functions [32]. The proposed DLCNN technique approximates the nonlinearity mapping that has a 
probabilistic network 7g (a|s) employed through a DLCNN with the parameter 0. The contender solution is 
sampled from the network probability, which means uw ~ 7g(als). The network’s loss function for 
continuous solution action domain. The symbols employed here are recognizable from traditional DLCNN 
learning for simplicity, and Figure 2 shows how they apply to a particular solution problem. 


Sampled solution {û, ~ z, (a | s, )} {û,} 
-Circuit condition in 
o (a |S ) Fá \ 
Network tt. Governing equation 
d Ho Oo ` eqnl(t) 
eqn2(t) 


Boundary conditon 
2 <t <25 and 0stS15 
Initial condition 
initial current I1(0) 2A 
r ` initial charge O(0) 2C 7 
ee TE i a e Se 


State (sampled computational domain Q) 1 5o 


Figure 2. DLCNN framework for the DEs of the adopted two-loop RLC circuit 


This research uses a DLCNN with six convolution layers and one fully linked layer. A max pooling 
layer, activation layer (linear rectified unit, or ReLU), and normalization batch layer are present in every 
convolution layer, with the exception of the final convolution level, when an average pooling level is utilized 
in place of the maximum layer. In the output layer, softmax activation is present. Figure 3 displays the 
screenshot for the specifics of the created network layers. The training configuration are: i) every nine 
epochs, the learning rate is decreased by a factor of 10; ii) setting the initial learning rate to 2e-2; iii) setting 
the maximum number of epochs to 12; iv) using a stochastic gradient descent with momentum (SGDM) 
solver with a mini-batch size of 256; and v) plotting the training progress. 
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trainedNet 23H OA 00 
Analysis date: 01-Nov-2022 14:58:02 s i gs S 
ANALYSIS RESULT @ 
Name Type Activations Learnables 
e Input Layer 
Input Layer Image Input 1x1024x2 - 
t 1x1024=2 images 
$ CNNI CNN1 Convolution 1x1024x16 Weights 1x8x2x16 
16 1*8*2 convolutions wit Bias 1x1x16 
t 
BN1 Batch Normalization |1x1024x16 Offset 1x1x16 
$BN1 Batch normalization with 1 Scale 1x1x16 
t $ MaxPool5 Boy ReLU 1x1024x16 - 
1 
grau t MaxPool1 Max Pooling 1x512x16 - 
t 1x2 max pooling with strid. 
© CNNG 
$ MaxPoolt CNN2 Convolution 1x512x24 Weights 1x8x16x24 
a t 24 1x8x16 convolutions wi Bias 1x1x24 
1 $ BNG BN2 Batch Normalization |1x512x24 Offset 1x1x24 
ocni li Batch normalization with 2 Scale 1x1x24 
t ReLU2 ReLU 1x512x24 - 
@ ReLU6 ReLU 
$ BN2 t MaxPool2 Max Pooling 1x256x24 Š 
1*2 max pooling with strid 
1 @AP1 
CNN3 Convolution 1x256x32 Weights 1x8x24x32 
$ ReLU2 tł 32 1*8=24 convolutions wi. Bias 1x1x32 
t e@FCi BN3 Batch Normalization | 1x256x32 Offset 1x1x32 
Batch normalization with 3 Scale 1x1x32 
© MaxPool2 a 
ReLU3 ReLU 1x256x32 
i $ SofiMax Hat 
$ CNN3 t MaxPool3 Max Pooling 1x128x32 - 
1x2 max pooling with strid. 
t 
© Output CNN4 Convolution 1x128x48 Weights 1x8x32x48 
@BN3 48 1*8=32 convolutions wi Bias 1x1x48 
t BN4 Batch Normalization | 1x128x48 Offset 1x1x48 
Batch normalization with 4 Scale 1x1x48 
TSAN ReLU4 ReLU 1x128x48 z 
t ReLU 
Š MaxPool3 MaxPool4 - Max Pooling 1x64x48 
1x2 max pooling with strid. 
t CNN5 Convolution 1x64x64 Weights 1x8x48x64 
64 1x8x48 convolutions w i 
saw 64 1x848 convolutions wi Bias 1x1x64 
BNS Batch Normalization | 1x64x64 Offset 1x1x64 
7 Batch normalization with 6 Scale 1x1x64 
qen ReLU5 ReLU 1x64x64 s 
t ReLU 
MaxPool5 Max Pooling 1x32x64 - 
$ ReLU4 1x2 max pooling with strid 
t CNN6 Convolution 1x32x96 Weights 1x8x64x96 
96 1=8x64 convolutions wi. Bias 1x1x96 
$ MaxPool4 
BN6 Batch Normalization | 1x32x96 Offset 1x1x96 
t Batch normalization with 9 Scale 11x96 
© CNN5 ReLU6 ReLU 1x32x96 - 
ReLU 
! AP1 Average Pooling 1x1x96 a 
$ BN5 1x32 average pooling with 
t FC1 Fully Connected 1x1x11 Weights 11x96 
11 fully connected layer Bias 11x1 
 ReLUs SoftMax Softmax 1x1x11 5 
softmax 
Output Classification Output |1x1x11 - 
crossentropyex with '16QA 


Figure 3. Screenshot for the details of the developed network layers 


3. RESULTS AND DISCUSSION 


3.1. Equations’ solving 


The Laplace transform is calculate for eqn1 and eqn2 as: 


eqniLT = laplace(eqni,t,s) = 


R2(Q(0)—s laplace(Q(t),t,s)) _ _ (RitRa)laplace(I1(t),t,s) 
L L 


laplace(I,(t),t,s) — 1,(0) + 


eqn2LT = laplace(eqn2,t,s) = 


E 
Rzlaplace(l4(t),t,s) -z7 laplace(Q(t),t,s) 
s laplace(Q(t), t, s) — Q (0) = #4 EAR so o AT 


By substituting Laplace laplace (Q(t), t,s) and (I1(t), t, s) by the variables Q_LT and I1_LT we get: 


(0) + R2(Q(0)-Qirs) _ _ ILLT(R1+R2) 


eqn1LT = hırs- h ; i 
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Cc 
ira = 
eqn2LT = Qiurs —Q(0) = ILLTR2 _ LT S244 0) 


R2+R3 c (R2+R3) 
Solving the equations for Q_LT and I1_LT we get: 
L14(0)—R2Q(0)+CRzS+Ls? I, (0)—RzS? Q(0)+C L R28? 14(0)+ 


C L R38? 14(0)+C L R3s11(0)+C (R1+R2+LS+L RzI4(0)+R1R29(0) 
(s? +1)(R1+R2+LS+C L Rs? +C L R3S?+C Ry Rost C Ry R3S+C R3 R35 ) 


I1 -LT = 


+R1R3Q(0)+R2R3Q(0)+L R25? I4 (0)+L Ros? Q(0)+L R38? Q(0)+R1R25s° Q(0)+ (10) 
R1R3S*% Q(0)+R1R3S7 Q(0)+L Ros Q(0)+L R3s Q (0) 
(s2 +1)(R1+R2+LS+C L R28? +C L R3s?+C Ry R2S+ C Ry R3S+ C R3 Rgs ) 


Calculate Q and 71 by calculating the inverse Laplace transform of Q_LT and I1_LT. Then 
simplifying the answer and substituting the symbolic variables with the circuit element values as listed in 
Table 1 to get so: and Qsor: 


Table 1. The symbolic variables with the circuit element values 


Description Value 
R1 42 
R2 22 
R3 30 
L 1.6H 
C “4 F 


Initial current 71 (0) 2A 


Initial charge Q (0) 2C 


. (Vi7eit 
742529 V1761 ECON 
aie at cosh EÀ 2 ( a 
40 14195421 
200 cos (t) , 405sin(t) 
[sol = 
8161 8161 8161 
aie serpy 1109425 VI76T sinn 5 ‘) 
~40 | 
17377e cosh( a0 + 30600897 
924 sin(t 1055 cos(t 
Qsol = azt Smt) __ 08 B cos) 


8161 8161 8161 


3.2. Laplace transform results 
Splitting the steady state and transient terms for Qsol and I1sol results in: 


I1steadystate = azo sin ©) 
8161 8161 
a gran 742529 Vre sinn( ES) 
1612 cosh( 40 )- 14195421 
I1steadystate = 


8161 


The drawing shown in Figure 4 represents the transient (on the left) and steady-state (on the right) 
behaviors for the current 71sol (a) and the charge Qsol (b) over two unlike periods given by 2 < t < 25 and 
0 < t < 15 is shown in, while Figure 5 shows the transient in (a) and steady state in (b) terms. 


Current Current Charge Charge 


1 Sgal Steady State 


0.5 Steady State 0.5 


(a) (b) 


Figure 4. The transient and steady state behavior for the charge: (a) Qsol and (b) the current J1sol over two 
unlike time periods given by 2 < t < 25 and 0 < t < 15 respectively 
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Likewise, splitting Qsol into steady state and transient terms expresses how symbolic computations 
assist to problem analysis, as given by: 


Qsteadystate Z (- 1055 cos (t) 924 sin ©) 


8161 8161 


40 30600897 


sit coat 1109425 i7eisinn( OE) 
1737 ~ 40| cos ( + 


Qtransient = 
8161 


Current terms Charge terms 


1sol 


| 
Q, 


0 5 10 15 Sa 5 10 15 
Time (sec) Time (sec) 


(a) (b) 


Figure 5. The transient and steady state terms of; (a) /1sol and (b) Qsol 


3.3. Comparison results 

Figure 6 shows the comparison results between the Laplace transform (LT) and the developed DLCNN 
of the transient in Figure 6(a), where in Figure 6(b) is steady state. For behaviors for the current //so/ and charge 
Qsol over two unlike time periods given by 2 < t < 25 (Figure 6(c)) and 0 < t < 15 (Figure 6(d)). The 
modeling DEs of the current J1sol and the capacitor charging Qso1 are derived according to the symbolic 
variables with the circuit element values of Table 1. The charge and current initially decline exponentially. On the 
other hand, they oscillate over a long period of time. The terms “steady state” and “transient” are used to obtain 
these characteristics. The qualities of the result can be examined with a symbolic result, but this is not possible with 
a numeric result. Figure 4 demonstrated that examining Qsol and /1sol visually shows that these equations consist 
of many terms. The plotting of the behaviors over [0 15] is performed to determine their contributions. The plots in 
Figure 5 show that Qsol has two steady-state terms and a transient, while Ilsol has a steady-state term and 
transient. From visual inspection, it is noted that Qsol and /1sol are with an exp function term, which is assumed 
to cause the decay in transient exponentially. 


T Charge {[——[T soluti 
= LT solution —LT solution 0.05 0.5 args LT solution 
2 — DLCNN solution 2 —*— DLCNN solution 04 —*— DLCNN solution 
Transient 
as 0 03 
S S02 Steady State 
o 
-0.05 V V 0.1 
Steady State 
= TL solution 7 
-0.1 f Transient PLONY solution od 
-0.5 - > d 
s a 
0 t 10 15 5 10 15 20 25 5 10 ‘ 15 ) 5 
(a) (c) (d) 


Figure 6. The comparison of transient and steady-state behavior for (a) the charge Qsol and (b) the current 
I1sol over two unlike time periods given by; (c)2 < t < 25 and (d) 0 < t < 15 


The Laplace transform and deep learning CNN solutions show excellent consistency for the steady 
solutions and there is a good agreement between the obtained transients. Future work can be recommended 
when comparing the solving of such symbolic differential equations of the 2-loop RLC AC circuit with deep 
learning convolutional neural network [33], [34]. 
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4. CONCLUSION 

This work solves the DEs of a two loops RLC circuit of an alternating voltage source by using 
Laplace transform. The capacitor charge and current initially decline exponentially and they oscillate over a 
long period. The qualities of the transient and steady-state results can be examined with a symbolic result, 
which is not possible with a numeric result. Three major points can be concluded: i) the approach presented 
has successfully modeled and solved the symbolic differential equations of the currents and the capacitor 
charging of the 2-loop RLC AC circuit; ii) the visualization of the obtained results of the capacitor charging 
and inductor current shows that these equations consist of many terms that are assumed to cause the decay in 
transient exponentially; and 111) this work presented the DLCNN-based solution and compared the results 
with its corresponding Laplace transform solution. 
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